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Abstract 



Scalar and vector bosons in the background of one- dimensional nonminimal vector linear plus 
inversely linear potentials are explored in a unified way in the context of the Duffin-Kemmer- 
Petiau theory. The problem is mapped into a Sturm-Liouville problem with an effective singular 
oscillator. With boundary conditions emerging from the problem, exact bound-state solutions 
in the spin-0 sector are found in closed form and it is shown that the spectrum exhibits degen- 
eracy. It is shown that, depending on the potential parameters, there may or may not exist 
bound-state solutions in the spin-1 sector. 

Key words: Dumn-Kemmer-Petiau theory, nonminimal coupling, Klein's paradox 



PACS Numbers: 03.65.Ge, 03.65.Pm 



1 Introduction 



The first-order Duffin-Kemmer-Petiau (DKP) formalism PQ-[1] describes spin-0 and spin-1 par- 
ticles and has been used to analyze relativistic interactions of spin-0 and spin-1 hadrons with 
nuclei as an alternative to their conventional second-order Klein-Gordon and Proca counter- 
parts. The DKP formalism enjoys a richness of couplings not capable of being expressed in the 
Klein-Gordon and Proca theories [2]- El- Although the formalisms are equivalent in the case 
of minimally coupled vector interactions [7J-[S]> the DKP formalism opens news horizons as 
far as it allows other kinds of couplings which are not possible in the Klein-Gordon and Proca 
theories. Nonminimal vector potentials, added by other kinds of Lorentz structures, have al- 
ready been used successfully in a phenomenological context for describing the scattering of 
mesons by nuclei [TU]- [TT] . Nonminimal vector coupling with a quadratic potential [IB], with a 
linear potential [IH], and mixed space and time components with a step potential [20]- [21] and 
a linear potential [22]- [23] have been explored in the literature. See also Ref. [22] for a com- 
prehensive list of references on other sorts of couplings and functional forms for the potential 
functions. In Refs. [22j-[23j it was shown that when the space component of the coupling is 
stronger than its time component the linear potential, a sort of vector DKP oscillator, can be 
used as a model for confining bosons. 

The Schrodinger equation with a quadratic plus inversely quadratic potential, known as 
singular oscillator, is an exactly solvable problem [24j-[28j which works for constructing solv- 
able models of N interacting bodies [29]- [30] as well as a basis for perturbative expansions 
and variational analyses for spiked harmonic oscillators [M]-[36j. Generalizations for finite- 
difference relativistic quantum mechanics [37] as well as for time- dependent parameters in the 
nonrelativistic version have also been considered [38j-[39j. 

The main purpose of the present article is to report on the properties of the DKP theory 
with nonminimal vector linear plus inversely linear potentials for spin-0 and spin-1 bosons in 
a unified way. The problem is mapped into an exactly solvable Sturm-Liouville problem of 
a Schrodinger-like equation. The effective potential resulting from the mapping has the form 
of the singular oscillator potential. The Schrodinger equation with quadratic plus inversely 
quadratic potential is indeed an exactly solvable problem. Nevertheless, only positive coeffi- 
cients are involved in the well-known solution. Because we need solutions involving a repulsive 
as well as an attractive inverse- square term in the effective potential, the calculation including 
this generalization with proper boundary conditions at the origin is presented. The closed form 
solution for the bound states is uniquely determined with boundary conditions which emerge 
as a direct consequence of the equation of motion and the normalization condition, and do not 
have to be imposed. It is shown that the spectrum exhibits degeneracy. It is also shown that 
the existence of bound-state solutions for vector bosons depends on a too restrictive condition 
on the potential parameters. 

It should be mentioned that a somewhat less general sort of problem, only with the space 
component of a nonminimal vector potential (erroneously called pseudoscalar potential), has 
already appeared in the literature [3D"] . 

2 Nonminimal vector couplings in the DKP equation 

The DKP equation for a free boson is given by [I] (with units in which h — c — 1) 

(i/3^-m)V = (1) 
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where the matrices (3^ satisfy the algebra (3 fl f3 u f3 x + (5 x f3 v f3^ = g fiu f3 x + g Xv (5^ and the metric 
tensor is g^ u = diag (1, —1, —1, —1). That algebra generates a set of 126 independent matrices 
whose irreducible representations are a trivial representation, a five-dimensional representa- 
tion describing the spin-0 particles and a ten-dimensional representation associated to spin-1 
particles. The second-order Klein-Gordon and Proca equations are obtained when one selects 
the spin-0 and spin-1 sectors of the DKP theory. A well-known conserved four-current is given 
by J M = ip(3 ll 'ip/2 where the adjoint spinor -0 is given by -0 = ip^rf with rf = 2(3° (3° — 1. The 
time component of this current is not positive definite but it may be interpreted as a charge 
density. Then the normalization condition J dr J° = ±1 can be expressed as 

J dr#V = ±2 (2) 

where the plus (minus) sign must be used for a positive (negative) charge. 

With the introduction of nonminimal vector interactions, the DKP equation can be written 

as 

(z/3^-m-z[P,/^)^ = (3) 

where P is a projection operator (P 2 = P and P^ = P) in such a way that ip[P, P^ip behaves 
like a vector under a Lorentz transformation as does ifj^ip. Once again d^J^ = [22]. If the 
potential is time-independent one can write ip(r,t) = 0(r*) exp(-iEt) , where E is the energy 
of the boson, in such a way that the time-independent DKP equation becomes 

[P°E + iP% - (m + i[P, p^AJ] = (4) 

In this case J M = 0/3^0/2 does not depend on time, so that the spinor describes a stationary 
state. 

The DKP equation is invariant under the parity operation, i.e. when r — > if A 

changes sign, whereas A Q remains the same. This is because the parity operator is V = 
exp(i5p)P rj°, where 5p is a constant phase and P changes r into — r . Because this unitary 
operator anticommutes with (3 l and [P, f3 l ], they change sign under a parity transformation, 
whereas 0° and [P, which commute with remain the same. Since 5p = or 5p = 
7r, the spinor components have definite parities. The charge-conjugation operation can be 
accomplished by the transformation ip — > tp c = Cip = CKip, where K denotes the complex 
conjugation and C is a unitary matrix such that C/3 M = —(5^C. Meanwhile C anticommutes 
with [P, /3 M ]. The matrix that satisfies these relations is C = exp (i5c) i] -)] 1 . The phase factor 
exp (iSc) is equal to ±1, thus E — > —E. Note also that — > — J 7 ', as should be expected for a 
charge current and the charge-conjugation operation entails no change on A^. The invariance 
of the nonminimal vector potential under charge conjugation means that it does not couple to 
the charge of the boson. In other words, A^ does not distinguish particles from antiparticles. 
Hence, whether one considers spin-0 or spin-1 bosons, this sort of interaction can not exhibit 
Klein's paradox. 

For the case of spin 0, we use the representation for the ^ matrices given by [H] 
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where 



P2 
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1 
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(6) 



0-1 




0, and are 2x3, 2x2 and 3x3 zero matrices, respectively, while the superscript T designates 
matrix transposition. Here the projection operator can be written as [5] P = (/^/^ — 1) /3 = 
diag (1, 0, 0, 0, 0). In this case P picks out the first component of the DKP spinor. The five- 
component spinor can be written as ip T = (?/>i, -0s) in such a way that the time-independent 
DKP equation for a boson constrained to move along the X-axis, restricting ourselves to 
potentials depending only on x, decomposes into 
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For the case of spin 1, the (3^ matrices are 
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-i£ijk, e% are the 1x3 matrices (ej 



\ x - = 5ij and 



where Sj are the 3x3 spin-1 matrices (sj) 
= (0 0), while I and designate the 3x3 unit and zero matrices, respectively. In this 
representation P = (3^(3^ — 2 = diag (1, 1, 1, 1, 0, 0, 0, 0, 0, 0), i.e. P projects out the four upper 
components of the DKP spinor. With the spinor written as ip T = (ipi, ...,if)i ), and partitioned 

as 

^8 = 
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the one- dimensional time-independent DKP equation can be expressed as 

4>\1 = ^ (£ + !<t-4o) # (12) 



where a is equal to + or — , and 



kl = E 2 -m 2 + Al \-A\ l + a^ (13) 



Now the components of the four-current are 

m t « m. * ' n.T. 



It is of interest to note that (p) (in the vector sector) obeys the same equation obeyed 
by the first component of the DKP spinor in the scalar sector and, taking account of (9) and 
(14), <pi and (pf are square-integrable functions. Given that the interaction potentials satisfy 
certain conditions, we have a well-defined Sturm-Liouville problem and hence a natural and 
definite method for determining the possible discrete or continuous eigenvalues of the system. 
We also note that there is only one independent component of the DKP spinor for the spin- 
sector instead of the three required for the spin-1 sector, and that the presence of a space 
component of the potential might compromise the existence of solutions for spin-1 bosons when 
compared to the solutions for spin-0 bosons with the very same potentials (Aq and A±). This 
might happen because the solution for this class of problem consists in searching for bounded 
solutions for two Schrodinger equations. It should not be forgotten, though, that the equations 
for and (f)\~^ are not independent because the energy of the boson, E, appears in both 
equations. Therefore, one has to search for bound-state solutions for and <p\ ' with a 
common energy. This amounts to say that the solutions for the spin-1 sector of the DKP 
theory can be obtained from a restricted class of solutions of the spin-0 sector. 



3 The linear plus inversely linear potential 

Now we are in a position to use the DKP equation with specific forms for vector interactions. 
Let us focus our attention on potentials in the linear plus inversely linear form, viz. 

A 2 \ \ , 90 A 2 i 9l I-, r\ 

A = m uj \x\ + — — , Ai = m uj\X H (15) 

\x\ X 

where the coupling constants, u , uji,g Q and g 1; are real dimensionless parameters. Our problem 
is to solve (J7J) and ffTiZj) for <p x (in the scalar sector) and (j)^ (in the vector sector), and to 
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and 



with 



determine the allowed energies. Although the absolute value of x in A is irrelevant in the 
effective equations for 4>i and <p^f\ it is there for ensuring the covariance of the DKP theory 
under parity. It follows that the DKP spinor will have a definite parity and will be a 
genuine four-vector. In this case the first equations of (J7J) and f|T2|) transmutes into 

H a $ v = e a $ a (16) 

where $ CT is equal to 4>% for the scalar sector, and to </>j for the vector sector, with 

1 d 2 

CT 2m dx 2 ° ^ 

E 2 -m 2 + m 2 [<jui + 2 {u g - ^i9i)] , . 

E ° = 2^ (18) 

1 n 

K = -mn 2 x 2 + ^ (19) 

2 x l 

il =m Iuj 1 - u Q ) , a a = ± (20) 

2m 

The set ffT6"j) - f[2T)j) . with a a = (a CT 7^ 0) and Q 2 > 0, is precisely the Schrodinger equation 
for the nonrelativistic nonsingular (singular) harmonic oscillator. For a a < and Q 2 = 
the effective potential has also a form that would make allowance for bound-state solutions 
with e a < 0. All the remaining possibilities for a a and Q 2 make the effective potential either 
everywhere repulsive or there appears an attractive region about x = between potential 
barriers (a a < and Q 2 < 0) which lead to boson tunneling. It is well known that the 
singularity at x = for a a < menaces the boson to collapse to the center [21] so that the 
condition a a > a CI ± t = — 1/ (8m) (with Q 2 > 0) is required for the formation of bound-state 
solutions. We shall limit ourselves to study the bound-state solutions. 

Since V a is invariant under reflection through the origin (x — > —x), eigenfunctions with 
well-defined parities can be built. Thus, one can concentrate attention on the positive half-line 
and impose boundary conditions on Q a at x = and x = 00. Continuous eigenfunctions on the 
whole line with well-defined parities can be constructed by taking symmetric and antisymmetric 
linear combinations of These new eigenfunctions possess the same eigenenergy, then, in 
principle, there is a two- fold degeneracy. As x — > 0, the solution behaves as Cx s<T , where C ^ 
is a constant and s CT is a solution of the algebraic equation 



s a (s a -l) -2ma a = (21) 

viz. 

s CT = - (1 ± VI + 8ma a ) (22) 

where s a is not necessarily a real quantity. Nevertheless, J 1 for a stationary state, as expressed 
by fl9]) and f TH|) . is the same at all points of the X-axis and vanishes for a bound-state solution 
(because vanishes as x — > 00), so we demand that s ff G M and so a a > ct crit . Normalizability 
of $0- also requires s a > —1/2 and due to the two-fold possibility of values of s a for a crlt < 
a a < 3/(8m), it seems that the solution of our problem can not be uniquely determined. This 
ambiguity can be overcome by a regularization of the potential. Following the steps of Ref. 
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we replace V a (x) by V a (x ) for x < x ~ and after using the continuity conditions for 
$ CT and d&ujdx in the cutoff we take the limit Xq — > 0. It turns out that the solution with the 
lesser value of s a is suppressed relative to that one involving the greater value as x — > 0. Thus, 
the minus sign in fl22|) must be ruled out for a a ^ in such a way that $ CT (0) = 0. Hence, 
s a = or s a > 1/2. The homogeneous Dirichlet boundary condition ($ CT (0) = 0) is essential 
whenever a a ^ 0, nevertheless it also develops for a a = when s a — 1 but not for s CT = 0. 
The continuity of § a at the origin excludes the possibility of an odd-parity eigenfunction for 
s = 0, and effects on the even-parity eigenfunction for s = 1 are due to the continuity of its 
first derivative. Effects due to the potential on d§ a /dx in the neighbourhood of the origin can 
be evaluated by integrating (Tl6|) from —5 to +5 and taking the limit 5 — > 0. The connection 
formula between d<& a /dx at the right and d^ a jdx at the left can be summarized as 



hm 



so that 



d$ r 

hm 



<5^o dx 



5->o dx 



x=+s 2mCa a 



-s 



•^s r+S 
lim / 

s^ols x 



2m« (T lirn^ / dx — (23) 



\x\* 



Sa — 1 



x=+S 

, for symmetric with a a ^ 

x=-8 ^ 

0, otherwise 



Therefore, the eigenf unctions have a first derivative continuous at the origin, except for s a < 1 
when the eigenfunction is not differentiable at the origin. Thus, the bound-state solutions for 
the singular oscillator potential are two-fold degenerate: for a same energy there is an odd- 
parity solution with a continuous first derivative, and an even-parity solution with a continuous 
first derivative if s a > 1 or a discontinuous first derivative if s a < 1. There is no odd-parity 
solution for s a = and there develops the homogeneous Neumann condition (cf$ a /dx\ x=0 = 0), 
and there is no even-parity solution for s a — 1. Thus, the bound-state solutions for the 
nonsingular oscillator potential are nondegenerate. 

We shall now consider separately two different possibilities for Q, namely Q = and Q ^ 0. 

3.1 Q = 

Defining 

z a = 2^—2me a x (25) 
for x > 0, and using the set f JT6|) - fr20|) . one obtains a special case of Whittaker's differential 



equation 

d 2 ^rr ( 1 2ma n 



- )$* = (26) 



dzl V 4 zl 

The normalizable asymptotic form of the solution as z a — >■ oo is e~ Za l 2 with z a > 0. Notice that 
this asymptotic behaviour rules out the possibility e a > 0, as has been pointed out already 
based on qualitative arguments. The exact solution can now be written as 

$ CT = z s /w(z a )e~ z ^ 2 (27) 

where w is a regular solution of the confluent hypergeometric equation (Kummer's equation) 

d 2 w .dw n , nnS 

Z ^~TT + (At _ z a)-j a a w = (28) 

dz£ dz a 
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with the definitions 



g<j s CT , b a 2s CT 



(29) 



The general solution of ( 1281) is expressed in terms of the confluent hypergeometric functions 
(Kummer's functions) \Fi(a a , b a , z a ) (or M(a a ,b a , z a )) and 2 -£o( a o-> 1 + a a ~ b a ,—l/za) (or 
U(a a ,b a ,z a )): 

w = A a i F\ ( , + B a z~ aa 2 Fo{a a , 1 + a a - b a , — -), b a ^ -n CT (30) 

where n a is a nonnegative integer. Due to the singularity of the second term at z a = 0, only 
choosing B a = gives a behavior at the origin which can lead to square-integrable solutions. 
Furthermore, the requirement of finiteness for & a at z a = oo implies that the remaining 
confluent hypergeometric function (iFi(a a ,b a , z a )) should be a polynomial. This is because 
iFi(a CT , b/j . z a ) goes as e Za as z a goes to infinity unless the series breaks off. This demands that 
a a = —n a , where n a is also a nonnegative integer. This requirement combined with f l2"§|) implies 
that the existence of bound-state solutions for pure inversely quadratic potentials (|a>i| = |u; |) 
is out of question. 



3.2 n^o 

As for Q 7^ 0, it is convenient to define the dimensionless quantity £, 

( = m\/ffi 2 (31) 
and using ( TT6|) - (|20|) . one obtains the complete form for Whittaker's equation 



^ + l^+/ r J^_f_^$ -o (32) 

The normalizable asymptotic form of the solution as £ — > oo is given by e~^ 2 only for f2 2 > 
(|cji| > | oj I ) - The solution for all £ can now be written as 

$ ct = £» °l 2 w{£)^' 2 (33) 
where u>(£) is a regular solution from 

*d 2 w ,, ^dw „ 

^ + ( 6 "-0-^--a^ = (34) 

with 

fl<T ~ 2 2|fi| 

(35) 

6 CT = s a + l/2 

Then, w is expressed as i-Fi(a CT , b a , £) and in order to furnish normalizable <& CT , the confluent 
hypergeometric function must be a polynomial. This demands that a a = —n a and b a ^ —n a . 
Note that i-Fi(— n CT , £) is proportional to the generalized Laguerre polynomial lit _1 '(0i a 
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polynomial of degree n a . The requirement a a = —n a combined with the top line of ( )35|) . also 
implies into quantized effective eigenvalues: 

2n a + s a + - ) \Q\, n a = 0,1,2,... (36) 



with 

$ CT oc ^ /2 e~ i/2 L^-W (0 (37) 



In summary, only for Q 2 > (e a > 0) and a a > « crit the potentials (fT5|) are able to furnish 
bound states. The spectrum is purely discrete and the solution is expressed in an exact closed 
form. It is instructive to note that s a = or s a = 1 for the case a a = and the associated La- 
guerre polynomial Ln~ 1/2) (0 and Ln^ 1/2) (0 are proportional to H 2na (a/T) or £,~ 1/2 H 2na +i {\fCj-, 
respectively [23]. Therefore, the solution for the effective nonsingular harmonic oscillator can 
be succinctly written in the customary form in terms of Hermite polynomials: 

n„ + ij|ft|, n a = 0,1,2,... (38) 



(v 7 ?) (39) 



with $ CT proportional to 

The preceding analyses shows that the effective eigenvalues for the scalar sector are equally 
spaced with a step given by 2|f2| when the effective potential is singular at the origin. It is 
remarkable that the level stepping is independent of the sign and intensity of the parameter 
responsible for the singularity of the potential. We also note that the effective spectrum varies 
continuously with a + as far as a + 7^ 0. When the effective potential becomes nonsingular 
(a+ = 0) the step switches abruptly to There is a clear phase transition when a + — > 
due the disappearance of the singularity. In the limit as a a — > the Neumann boundary 
condition, in addition to the Dirichlet boundary condition always present for a + 7^ 0, comes 
to the scene. This occurrence permits the appearance of even Hermite polynomials and their 
related eigenvalues, which intercalate among the pre-existent eigenvalues related to odd Her- 
mite polynomials. The appearance of even Hermite polynomials makes $+(0) 7^ and this 
boundary condition is never permitted when the singular potential is present, even though a + 
can be small. You might also understand the lack of such a smooth transition by starting 
from a nonsingular potential (aq_ = 0), when the solution of the problem involves even and 
odd Hermite polynomials, and then turning on the singular potential as a perturbation of the 
nonsingular potential. Now the "perturbative singular potential" by nature demands, if is 
either attractive or repulsive, that $+(0) = so that it naturally kills the solution involving 
even Hermite polynomials. Furthermore, there is no degeneracy in the spectrum for a + = 0. 

Now we move on to match a common energy to the spin-1 boson problem. The compatibility 
of the solutions for $ + = cj)^ and $_ = (j)^ demands that the quantum numbers n + and n_ 
must satisfy the relation 

n + — n_ = R(u , lux) (40) 

for the nonsingular potential, and 

1 



n+ - n_ = - 



R(u Q , wi) - a/1 + 8ma + + a/1 + 8ma_ (41) 



S 



for the singular potential. Here 



and sgn(cJi) stands for the sign function. For the nonsingular potential \n + — n_| varies from 
1 to infinity (as |wo/^i| varies from to 1), whereas for the singular potential \n + — n_| varies 
from the minor to the major value between \R(uo,ui)\ /2 and \R(ujq,uj\) ± \/&gi \ /2 (the ± 

sign corresponding to a± = a crit ). These constraints on the nodal structure of <fi^ and <f)\ ^ 
dictate that acceptable solutions only occur for a restricted number of possibilities for the 
potential parameters and that no solution should be expected for gi < 0. More than this, 
these constraints exclude a few low-lying quantum numbers for a given set of possible solutions 
for the simple reason that n + and n_ are nonnegative integers. In conclusion, the possible 
bound-state solutions for spin-1 bosons exist only for a finite set of potential parameters. 

To have a better understanding of the spectrum, we plot in Figured] the effective spectrum 
for the four low-lying levels as a function of a + , for fixed Uq and go . From this figure we see the 
phase transition for the spin-0 spectrum at a + = and the severe restriction on the solutions 
for the spin-1 spectrum. Notice that no bound-state solution for the spin-1 spectrum exists if 
a + exceeds \R(uq,ui)\/2 and that not all the bound states have a nodeless <3>+ for the ground 
state. 



4 Conclusions 

We approached the problem of a particle in a (3+l)-dimensional world despite the restriction 
to the one-dimensional movement. The motion on axis allow us to explore the physical con- 
sequences of the nonminimal vector coupling in a mathematically simpler and more physically 
transparent way. There is no angular momentum so that there is no spin-orbit coupling. The 
more general potential matrix ensuring a conserved four-current can be written in terms of 
well-defined Lorentz structures. For the spin-0 sector there are two scalar, two vector and two 
tensor terms [5], whereas for the spin-1 sector there are two scalar, two vector, a pseudoscalar, 
two pseudovector and eight tensor terms [6J. There is no pseudoscalar potential in the spin-0 
sector of the DKP theory. In fact, the space component of a nonminimal vector potential is 
used in Ref. [40]. We described spin-0 particles by a five-component spinor, of which only one 
component is independent instead the two required in Ref. [ID]. Also, the number of degrees of 
freedom described by the ten-component spinor in the spin-1 sector reduces correctly to three. 
By considering space plus time components of a nonminimal vector potential, we showed that 
scalar and vector bosons in the background of one-dimensional linear plus inversely linear po- 
tential is mapped into a Sturm-Liouville problem with an effective singular oscillator. The 
existence of solutions for spin-0 bosons requires \cui\ > \oj \ and g\ (g\ + a) — g\ > —1/4. In 
words, the time component of the nonminimal vector potential (Aq) alone can not hold bound 
states, and the presence of the time component in the inversely linear potential lessens the 
interval of parameters which provides confinement. As for spin-1 bosons, confining solutions 
requires g\ > and when compared to the solutions for spin-0 bosons with the very same 
potentials (Aq and A{) they can only be obtained from a restricted class of solutions of the 
spin-0 sector. It may be salutary to note that, because the differential equation ffTBj) has a 
singularity at the origin for a a ^ 0, one could expect the existence of singular solutions for <3> CT 
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so that we are not free to suppose that s a > 0. Due to the fact that whether they are even 
or odd, disappears at the origin for a a 7^ this one- dimensional quantum-mechanical problem 
is degenerate and H a in (TT7j) is a self-adjoint operator. 
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Figure 1: Effective spectrum (e = e + /m and n = n + ) for spin-0 bosons as a function of 
a = a + 7^ (full line) with uji — m — 1, ujq — 1/2 and g = 0. The circles and boxes stand 
for even and odd solutions for the case a + = 0, respectively. The effective spectrum for spin-1 
bosons are represented by the diamond symbols for n + — n_ = +2, +1, 0, —1 corresponding to 
a = -0.115, -0.047, 0.083, 0.275 respectively. 
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